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Abstract 

We study the dynamics of two-dimensional coherent structures in planetary atmospheres and oceans. 
We derive the Zakharov-Kuznetsov equation for large scale motion from the barotropic quasigeostrophic 
equation in a weakly nonlinear, long wave approximation. We consider coherent structures emerging 
out of an instability caused by a narrow jet-like meanflow. We use multiple scale analysis combined 
with asymptotic matching. 
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1 Introduction 

Planetary atmospheres and oceans are strongly turbulent media. However, highly ordered coherent struc- 
tures arise in a process of self-organization, and dominate the dynamics on slow temporal and large spatial 
scales. 

The spontaneous appearance of coherent structures is a characteristic of two-dimensional fluid flows. 
The basic underlying structure of these flows is linked to the existence of two quadratic, positive definite 
invariants, energy and enstrophy. In spectra of two-dimensional turbulence one observes two different 
cascades associated with these conserved quantities; a direct enstrophy cascade towards small spatial 
scales, and an indirect energy cascade towards larger spatial scales. It is the latter which gives rise to 
vortex merging leading to larger and larger vortices. In this Letter we view the problem as one of weakly 
nonlinear hydrodynamic stability rather than turbulence phenomenology. 

Most vortices are monopolar, but dipole and even tripole vortices can also appear spontaneously. Monopole 
vortices are mainly created by shear flow instabilities whereas dipole vortices typically appear when some 



1 



additional forcing is applied to the flow. 

The richness and complexity of two-dimensional flows and the simultaneous presence of motion on very 
different temporal and spatial scales makes a direct analysis of the basic equations of motion very 
difficult. The presence of rotation reinforces the two-dimensional character in accordance with the 
Taylor- Proudman theorem, but rotation can also introduce baroclinic instability. The latter is a three- 
dimensional feature and, thus supports a direct energy cascade towards small scales. Hence, the dynamics 
is determined by competing two-dimensional and three-dimensional processes |^nEll^- To study vor- 
tices in geophysical fluid dynamics the primitive equations are further reduced by approximations which 
allow to focus on temporal and spatial length scales of vortices EH] ■ If additionally baroclinic pro- 
cesses are excluded a further simplification can be made. The dynamically important variable is the so 
called potential vorticity q. The resulting quasigeostrophic barotropic vorticity equation 

D + f{y) 

—q = where q = — (1) 

describes large scale motion on a slow time scale. Here ip is the stream function, f{y) describes the 
ambient rotation of the planet and H is the fluid depth. This equation was first derived by Charney 
and then independently by Obukhov [2R]| . In the context of low- frequency drift waves in magne- 
tized plasmas Equation JQ) is known as the Hasegawa-Mima equation ,11 . This equation has been the 
mathematical starting point for much of the research done on coherent structures and vortices. It sup- 
ports so called modons which arc localized soliton-like coherent solutions. Exact modon solutions were 
obtained by Larichev & Reznik [ISj for a stationary double-vortex solution which is antisymmetric in 
longitude. Extensions to more general solutions have been made and the spherical geometry of 

planets has been incorporated 1811 1^ . However, modons have the drawback that the potential vor- 
ticity is not a smooth function of the stream function, but may be multivalued. Therefore interest has 
grown in low-dimensional models, although a rigorous proof of existence of a low-dimensional attractor 
in quasigeostrophic systems is still an unsolved problem. Strictly speaking, one can only define a "slowest 
invariant manifold" since the small-scale events, i.e. the high-frequency and high-wavenumber pro- 
cesses, enlarge the Hausdorff dimension for the attractor without any convergence J^^. Nevertheless, in 
order to understand better the particular mechanisms involved in the formation and dynamics of vortices 
in geophysical fluid dynamics, it is useful to perform asymptotic techniques to derive reduced amplitude 
equations of the basic quasigeostrophic equations in a multiple scale analysis and study the derived model 
evolution-equations. The basic idea is that coherent vortices may be identified with solitary wave solu- 
tions of generic nonlinear dispersive wave equations. 

Most research has been done in the frame work of the Korteweg-de Vries equation |21[ I33L 1271 IIUL 1191 
I2(J[ I23L E] or in the framework of the Boussinesq equation jl2[ [T^ . While these models were helpful 
in describing and identifying mechanisms for atmospheric blocking, cyclogenesis, meandering of oceanic 
streams and the persistence of the Great Red Spot in the Jovian atmosphere, they are all one-dimensional 
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models with their obvious Umitations. 

In this Letter we will extend weakly nonlinear, long wave multiple scale analysis to two dimensions and 
derive the Zakharov-Kuznetsov equation 



The Zakharov-Kuznetsov (ZK) equation is one of two well-studied canonical two-dimensional exten- 
sions of the Korteweg-de Vries equation 17 ; the other being the Kadomtsev-Petviashvilli (KP) equation 
In contrast to the KP-equation, the ZK-equation has so far never been derived in a geophysical 
fluid dynamics context. For a derivation of the KP equations for internal waves, see Whereas the 
KP-equation is valid in isotropic situations, the ZK-equation is valid in anisotropic settings which is 
exactly the case for rotating fluids where the differential longitudinal dependence of the rotation rate 
causes anisotropy between the meridional and the longitudinal directions. Moreover, in contrast to the 
KP-equation the ZK-equation supports stable lump solitary waves. This makes the ZK-equation a very 
attractive model equation for the study of vortices in geophsyical flows. 

The Letter is organized as follows. In Section 2 we set up the barotropic vorticity equation and the mean 
flow configurations under consideration. In the beginning of Section 3 we will give a simple heuristic 
scaling argument based on the linearized barotropic vorticity equation to motivate why the ZK-equation 
is the generic two-dimensional nonlinear wave equation. In the remainder of Section 3 we will derive the 
ZK-equation in an asymptotic multiple scale analysis. Section 4 concludes the Letter with a discussion 
and an outlook on further research. 

2 Barotropic Quasigeostrophic Equation 

We shall use a non-dimensional coordinate system, based on a typical horizontal length scale Lq, a typical 
vertical scale Hq, and typical Coriolis parameter /o. A typical velocity U is taken to be the maximum of 
the mean current velocity and the time scale is given by U / Lq. If we separate the meridional meanflow 
U from the perturbation pressure fields p and use the Boussinesq approximation, we obtain the following 
equation for the non-dimensional perturbation pressure field j28j 



At + AAx - fiAAx - ^Axxx - (Axyy = . 




(2) 



where 



q = 



VV - Ftp , 



y 



13-Uyy+FU , 
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with Froudc number F and the Jacobian defined by J (a, b) = a^by — ayb^- We investigate a channel flow 
with a storm track superimposed on a constant meanflow confined at y — ±L (see Fig.l). The storm 
tracks may have a critical layer where U{y) = 0. Important is, as we will see, the non-vanishing slope at 
at least one boundary of the localized storm track. The boundary conditions are ip — const at y — ±00, 
and we require that the jet forms a transport barrier to the flow. 



-I ^ 
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y 
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Figure 1: Sketch of a typical mean flow. 



3 Nonlinear Wave Equation 
3.1 Linear Dispersion Relation 

Before we consider the weakly nonlinear, long wave approximation, we motivate our approach by looking 
at the linearization of equation ^ which yields 

l + U—]q+^^ = (3) 
dt dx J dx dy 

In terms of "0 = aoexpi(k{x — ct) + ly) we obtain the dispersion relation 

c = U ^ 

provided that the meanflow U{y) is constant. If wc focus on dynamics on the spatial and temporal time 
scales X — ex,Y = ey, i.e. small k, I, and T = e^t this is suggestive of the coupled Zakharov-Kuznetsov 
equation 

At + AAx - ^JiAAx - £.Axxx - CAxyy = . (4) 
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The reason why we expect an equation of the Zakharov-Kuznetsov type instead of the usual Kadomtsev- 
PetviashviUi type mostly encountered in fluid systems is the anisotropic character of caused by the 
/3-effect. 

3.2 Weakly Nonlinear Model 

We consider weakly nonlinear waves riding on a background meanflow. The meanflow consists of a 
constant part U„i and a strong but narrow Jetstream (see Fig. 1). The narrow storm track is located on 
a short meridional scale y. In the outer region the problem ||2Jl can be reduced to the linear problem © 
with constant meanflow In the interior the structure of the storm track does not allow for sinusoidal 
wave solutions but instead we will derive a nonlinear wave equation. In order for the nonlinear wave 
equation which is valid only in the inner region where the storm tracks are nonuniform, the inner solution 
has to be matched to the outer sinusoidal solution. 

3.2.1 Outer solution 

In the outer region where the meanflow is uniform and constant, ((SJ reduces to the simple linear equation 
^ for the streamfunction with constant coefficients j28| . 
There the solution of the streamfunction can be written as 

= a{t, T, X) sin(iy) + h{t, T, X) cos{lY) , (5) 

where / is the meridional wave number and is determined by the dispersion relation of the linearized 
model (jSjl. 

3.2.2 Inner solution 

In the interior of the storm track on the small scale y, the meanflow is not constant. We shall study 
weakly nonlinear long waves. We introduce the following scales, 

X = ex, Y = ey, T = e^t, 

^(X,y,T,y) = e2^(") + £3^(1) + eV^) ^ ... ^ 

where e is a small parameter, the inverse of which measures the large scales of the disturbance. Next, 
we rescale the parameters F e^F and /3 The scaling of the Froudc numbers implies that our 

model is valid for situations where the internal Rossby radius of deformation is of the order of the long 
horizontal scale. Further, the scaling of f3 implies that Qy « ^Uyy at the lowest order. The boundary 
conditions we use are ip — constant ai y — ±oo and Uyipx — Uipxy at y = ±L which simply states that 
there is no transport of fluid across the jet-stream. 
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Substituting this scaling into equation Q yields, 

= {e'^dr + eUdx){e^dxxi'^"'> + e^dyyi/°'> + 2e^dyY^^°'> + e^dyyi^^^^ " e^F^^°^ 
+ e^'dyy^'^^l + e^dyyip'^^^ + 2e%Yi^^'^) 
+ e^^f{e^(3-Uyy + e^FU) + eV^'^(-f/,,) + e'^kfi-Uyy) 

Counting the orders of e we obtain to the lowest order, 0{£^) 

JO^'-x^ = , (6) 

where 

^ = ^ [Udy - Uy] . (7) 
We look for an amplitude equation, ie we want to write 

i;^°\X,Y,T,y) ^ A{X,Y,T)^^°\y) (8) 
and seek an evolution equation for the slowly varying amplitude A{X, Y, T). We easily find 

where aa is a constant of integration. If we allow for zero meanflow within the narrow jet region we need 
to impose ao = 0. We summarize the solution of equation © 

^^"Hx,Y,T,y) = A{X,Y,T)U{y). (10) 

The meridional structure on the small scale y of ■0 is entirely determined by the mean currents at the 
leading order. 

At the next order, 0{e^), we obtain a linear inhomogeneous equation for ip'^^K 



which can be written, using as 



•Jx = ~^^Vy ^XY ■ 

This equation is again solved by the method of variation of parameters and we obtain 
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with 



where 



-L U^y' 

We note that the higher order term ■0'^^^ is slaved to the ip'"^^ term and the dynamics of the corresponding 
amphtude equation which wiU be derived shortly. For the same reasons as above we set ai = and 
obtain 

V'^^^ ^U{l-y - L)Ay . 
The O(e^) terms give us an evolution equation for the amplitude A. We obtain 

U^Xyy-Uyy^X+^U^yYX 

+ - €'^'xly + + FU^'x'} - , 

which, using JHl, can be written as 

C^'^P = -G , (11) 

G = ip'^jAT + U^^''HAxxx+AxYY)-FU^^"'>Ax 
+ ^(°) (/3 + FU) Ax + (^(°V^"J, " ^P^g^) AAx 

+ 2Uv\^^AxYY. (12) 

To assure boundedness of the solutions of (|ll|l we have to require a solvability condition in form of a 
Fredholm alternative. 

The homogeneous adjoint problem to equation pi(l may be written as 

£^ = 0, (13) 

with 

= 2Uydy + Udyy , 

where we have used the boundary conditions Uyipx = Uipxy at y = ±L. The adjoint eigenvalue problem 
\i'6\ has one trivial constant kernel mode (pi = const and one nontrivial, namely 

The nontrivial kernel mode (p2 has to be discarded because it does not satisfy the boundary condition. 
To see this note that the Fredholm alternative for the elliptic operator [7| together with the boundary 
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condition Uyipx = UipXy is equivalent to one for the operator £' — dy[Udyip] ~ Uydyip ~ Uyydyip with 
the boundary condition dyijj = at y = ±L. Also note that a non-zero kernel mode (f)i is only consistent 
with the boundary conditions if U{±L) ^ 0. 

The solvability condition is thus given by the trivial constant kernel mode 

Gdy = 0. (15) 



On substituting the expressions (jHJ with ip^^\y) — U{y) and l(T^ we obtain the desired amplitude 
equation for A, 



At + AAx - fJ-AAx - iAxxx - C,Axyy = , (16) 



where 



7C = {U\\-y-L)\\ , 

= -[U'ytL + lUUyyt^ , 
lA = - J PU dy . (17) 

We note that due to the last term of (|12|l the Zakharov-Kuznetsov equation is inhomogeneous in the 
sense ^ C- It is pertinent to mention that a nonzero C. requires a nonzero mean flow at at least one 
of the boundaries of the storm track. The coefRcients of the nonlinear terms fi require a non-vanishing 
slope at at least one boundary. The slope Uy and also Uyy at the boundaries of the jet y ~ ±L may be 
determined from a given meanflow configuration by averaging over a very short region, say y/e, where a 
sudden change of the constant mean flow Um to the jet occurs. 



3.3 Asymptotic Matching 

At the lowest order the inner solution © with ao = and the outer solution © have to be matched. The 
outer solution has been derived on the large scale Y whereas the inner solution and its associated ampli- 
tude equation, the Zakharov-Kuznetsov equation l|16(l . were derived on the short scale y. Henceforth we 
need to require that the asymptotic limit of the outer solutions for y — > coincides with the asymptotic 
limit of the inner solution for y ^ oo. The limit of the inner solution is = A{X, Y, T)Um- The limit 
of the outer solution is = b{X, Y, T). Hence we find h{X, Y, T) = A{X, Y, T)[/„„ which extends the 

dynamics of the Zakharov-Kuznetsov equations to the outer region. 
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4 Discussion 



We have derived the nonhnear dispersive Zakharov-Kuznetsov equation from the quasigeostrophic barotropic 
vorticity equation. It is weU known that the ZK-equation, ahhough it is not integrable by means of the 
inverse scattering transform, supports a family of steady-shape stable lump solitary waves, moving at 
an arbitrary velocity |85l I15| . These may help to describe two-dimensional coherent structures such as 
atmospheric blocking events, long lived eddies in the ocean or coherent structures in the Jovian atmo- 
sphere such as the Great Red Spot. The model is from an analytical point of view easier to treat than the 
full barotropic quasigeostrophic equation and its solutions do not exhibit multivalued potential vorticity- 
stream function relationships as modons do. 

Geophysical flow on large scales is widely accepted to be conservative. This allows for Hamiltonian de- 
scriptions of the flow on large scales. Our model also exhibits a Hamiltonian structure. Note that the 
momentum 

/oo 
dX dY 
-OO 

and the Hamiltonian with the Hamiltonian density 
are conserved. 

We have assumed a meridional meanflow U which consists of a constant part Um and a narrow localized 
storm track. Note that the jet stream may also be a narrow interface between two regions of meanflow 
with opposite flow direction. Such persistent shear layers exist between the zones and belts in the Jovian 
atmosphere. 

Analysis of the solutions of is planned. Their stability has to be numerically tested within the 
Zakharov-Kuznetsov system. The ZK-equation has been derived using asymptotic techniques and is as 
such an asymptotic limit to the barotropic quasigeostrophic vorticity equation. However, it is not clear 
that the same is true for the solutions. The solutions of the Zakharov-Kuznetsov equations do not neces- 
sarily have to be asymptotically close to the solutions of the full quasigeostrophic system, this is due to 
the lack of a centre manifold as discussed in the introduction. In further work we will test the approxi- 
mation of the solution numerically by taking solutions of the ZK-equation and testing their dynamics in 
the full quasigeostrophic system. 
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